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ABSTRACT 

We present the results of a deep, near-infrared, narrow band imaging survey at a cen- 
tral wavelength of 1.062 (FWHM=0.01^m) in the GOODS-South field using the ESQ 
VLT instrument, HAWK-I. The data are used to carry out the highest redshift search for 
[OII]A3727A emission line galaxies to date. The images reach an emission line flux limit (5(7) 
of 1.5 X 10^'^erg cm^^ s^\ additionally making the survey the deepest of its kind at high 
redshift. In this paper we identify a sample of [OII]A3727A emission line objects at redshift 
z~ 1 . 85 in a co-moving volume of 4 1 00 Mpc^ . Objects are selected using an observed equiv- 
alent width (EWoii) threshold of EW„/„ > 50A. The sample is used to derive the space density 
and constrain the luminosity function of [Oil] emitters at z=1.85. We find that the space den- 
sity (p) of objects with observed [Oil] luminosities in the range log(L[oii]) > 41.74 erg s^' 

is log(p)=— 2.45±0.14Mpc^^, a factor of 2 greater than the observed space density of [Oil] 
emitters reported at z^l.4. After accounting for completeness and assuming an internal ex- 
tinction correction of A//a=l mag (equivalent to A[o//] = 1.87), we report a star formation rate 
density of ^0.38±0.06 M© yr^' Mpc We independently derive the dust extinction of 
the sample using 24/im fluxes and find a mean extinction of Aro//]=0.98±0.11 magnitudes 
(A//a=0.52). This is significantly lower than the A//a=l (A[o//]=1.86) mag value widely used 
in the literature. Finally we incorporate this improved extinction correction into the star for- 
mation rate density measurement and report p* ^ 0.24 ± 0.06 M© yr^' Mpc 

Key words: galaxies:high-redshift, galaxies:luminosity function, galaxies:star formation, 
galaxies:distances and redshifts 



1 INTRODUCTION 

The volume-averaged star formation history is a fundamental prop- 
erty of the universe: its reliable determination will provide a pow- 
erful probe with which to explore the physics of galaxy formation 
and evolution. 

Much effort has been focused in this field and a variety of 
star formation rate (SFR) indicators at different redshifts have con- 
tributed to the cuiTent picture whereby star formation starts be- 
tween z~20-8, peaks at z~l-3 and then declines towards z=0. Al- 
though constrained within 30-50% below z-^l, the star formation 
rate is less well determined at higher redshifts, being uncertain up 
to a factor of 3 between 1 < z < 6 (Hopkins & Beacom, 2006). 

Variation between different measurements is primarily due to 
differences in sample selection, biases between different SFR in- 
dicators and underlying cosmic variance. Whilst using a combina- 
tion of different indicators has provided a qualitative description 
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of the evolution of the SFR, such biases make it diffi cult to prop- 
erly quantify the evolution. Indeed, as pointed out by iGeach et al.l 
(2008) , piecing together measurements from different indicators is 
no longer improving our understanding. 

To make progress, homogeneous indicators are needed, visible 
over wide redshift ranges. Although Ha remains the SFR indicator 
of choice, the [Oil] A 3727 doublet has a particular advantage over 
Ha in being visible to z~ 5 in the near-infrared, compared to the 
z~ 2.5 limit of Ha surveys. 

Up until recently, [Oil] has for the most part remained on the 
periphery of efforts to measure the SFR history of the universe. As 
a collisionally excited forbidden line, the [Oil] doublet is not di- 
rectly coupled to the UV ionising radiation and as such, the [OII]- 
SFR calibration is subject to scatter due to excitation variations re- 
lated to metal abundances and ionisation state (see the review by 
Kennicutt, 1998). Despite these limitations, the intrinsic [OII]/Ha 
variation is typically a factor of two or less over wide ranges of 
galaxy environments and abundances and [Oil] has been employed 
in a range of previous studies (see Hippelein et al. 2003; Hopkins 
2004; Kewley, Geller & lansen 2004; Ly et al. 2007; Takahashi et 
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al. 2007; Zhu, Moustakas & Blanton 2009) acting as a useful index, 
particularly for large statistical samples. 

Recent advances in the calibration of the [Oil] A 3727 dou- 
blet (see Kewley et al. 2004; Moustakas, Kennicutt & Tremonti 
2006; Kennicutt et al. 2009) have made it possible to use [Oil] with 
greatly improved precision, approaching that of more traditional 
SFR indicators. Addressing concerns over metallicity and excita- 
tion, Moustakas et al. (2006) found that the majority of variation 
between [Oil] and SFR is due to dust-reddening (derived from the 
Ha/Hj3 decrement) and that variations in metallicity and excitation 
are in fact second-order effects in most galaxies. 

Taking this further, Kennicutt et al. (2009) developed empiri- 
cal calibrations between [Oil] luminosity and SFR using weighted 
combinations of either Total Infrared (TIR), 24/xm or 8/im flux to 
correct for dust extinction. They report that for z=0 galaxies, the 
dispersion of [Oil] flux, corrected using their empirical relations, 
is equivalent to that of their corrected Ha samples, thus facilitat- 
ing the first reliable [OII]-derived measurements of the SFR to be 
made. 

In this paper, we apply these advancements in the calibration 
of [Oil] to the highest redshift narrow band survey for [Oil] emit- 
ters to date, concentrating on objects in the GOODS field at z=1.85. 
Using the new HAWK-I instrument on the ESO VLT facility, the 
survey covers a co-moving volume of ~ 4100Mpc-' to a depth of 
1.5 X 10^'^erg cm^^ s^'. 

In Section 2 we describe the data set, data reduction and 
cataloguing techniques used in the study. Section 3 describes the 
method for selecting emission line galaxies (ELGs) and how we 
remove emitters other than [Oil] from the sample. After taking 
into account the completeness, we compute the luminosity func- 
tion in Section 4. Section 5 looks at the number density evolu- 
tion of [Oil] emitters between z=0.8 and z=1.85. In Section 6, we 
convert the integrated [Oil] luminosity into a star formation rate, 
firstly using a standard extinction correction and secondly using 
the improved locally derived calibrations of Kennicutt et al. (2009). 
Throughout this paper, magnitude measurements are on the AB 
scale (m^5 = 48.60 — 2.51ogi()flux). A standard cosmology is as- 
sumed with SIm = = 0.7 and h = 0.70. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Observations 

This study utilises deep NIR data in two overlapping bands, a broad 
Y filter centred on 1.021/im and a narrow band filter at 1.060/im. 
Filter centres and FWHM values were calculated using the pre- 
scription in Pascual, Gallego & Zamorano (2007). Filter transmis- 
sion profiles are plotted in Figure [T] The data were obtained us- 
ing HAWK-I (Kissler-Patig et al. 2008), a NIR (0.85-2.5/xm) wide 
field, cryogenic imager on UT4 at the ESO VLT facility in Paranal, 
Chile. The HAWK-I focal plane is made up of 4 square Hawaii 
detectors separated by a cross-shaped gap of 15". Each filter has a 
2048 pixel width and a pixel scale of 0. 106"/pix. The HAWK-I field 
of view of is 7.5 x 7.5 arcmin. The data were collected in the Sci- 
ence Verification phase as part of program 60.A-9284(B) Fontana 
et al.: A deep infrared view on galaxies in the early Universe. A 
single pointing was taken in a region of the GOODS field centred 
on co-ordinates 3h 32m 29.0s, -27°44'28". 

Table[T]gives basic filter information, the average seeing of the 
images and the total exposure time in each filter. 
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Figure 1. Transmission curves of the Y and NB1060 titters used in this 
study. Filter centres and FWHM values are indicated by dashed and solid 
lines for the Y and NB1060 filters respectively. 



Table 1. Filter Statistics and Image Seeing. 



Filter 


Xc (Mm) 


AA (^im) 


Seeing (") 


Exposure Time 










(hr:min:s) 


NB1060 


1.0619 


0.0104 


0.74 


8:20:00 


Y 


1.0193 


0.1018 


0.57 


1:10:30 



2.2 Reduction and calibration 

The observations are made up of a series of spatially offset 
(dithered) exposures, comprising 141 x 30s exposures in Y and 
100 X 300s exposures in NB1060. The data was reduced using a 
pipeline specially developed for HAWK-I at the Cambridge As- 
tronomy Survey Unit that incorporates components of the VISTA 
Data Flow System (Irwin et al., 2004). The pipeline can be sum- 
marised as follows: 

Firstly, each exposure has its dark current and flat field instru- 
mental signature removed. The dark current was subtracted using a 
master dark image created from dark frames with the same expo- 
sure time as the image (about 100 dark frames were available per 
exposure time). 

The images were then flat fielded to remove the pixel to pixel 
quantum efficiency variation as well as the large scale vignetting 
profile. A master flat field for each filter was formed from exposures 
of the twilight sky, which were scaled to bring them to a common 
median background and then combined using a mean combina- 
tion algorithm with sigma-clipping (5cT). The telescope was moved 
slightly between each twilight flat exposure so that when twilight 
flat field images are combined, any remaining astronomical objects 
are removed by the rejection algorithm. 

Any pedestal scale factors between the individual detectors 
(due to differences in gain or average QE) were removed by nor- 
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Table 2. Source Extractor parameters - NB1060 source catalogue. 



Parameter 


Unit 


Value_^ 


Value_^ 


DETECT JVIINAREA 


pix 


8 




DETECT.THRESH 


a 


L5 




ANALYSIS.THRESH 


c 


L5 




BACK.SIZE 


pix 


64 


15 


BACKJILTERSIZE 


pix 


5 


3 


DEBLENDJSfTHRESH 


ADU 


32 




CLEAN 




Y 




CLEAN J-ARAM 




2.0 




PIXEL.SCALE 


"/pix 


0.1064 





" Parameters used for detectors 1, 3 and 4 

^ Changes made to standard parameter input for detector 2 to compensate for 
low-level background variations 



malising the images by the ensemble median of the background 
flux. 

Once the jittered exposures are dark and flat corrected, they 
are median combined to form a sky background image. The back- 
ground images are normahsed to a zero median so that the median 
background level in the corrected images is preserved. The back- 
ground corrected images from each jitter sequence are registered 
internally using visible sources and then shifted and combined. 

A world coordinate system is added by fitting sources detected 
on the image to the 2MASS point-source catalogue. Finally, all the 
stacks are combined together to form a single deep stack for each 
filter. 

The zero point of the Y images was determined by calibrating 
stellar counts onto ISAAC J and H images (Retzlaff et al. 2010), 
using synthetic J — H and Y — J colours generated bv lHewett et al.l 
( l2006h as reference. Y — NB1060 colours generated by the same 
synthetic codes were used to calibrate the NB 1060 images onto the 
Y band. 

The images were then scaled to a common zero point of 30.0 
(AB), accurate to 0.1 magnitudes. 



2.3 Cataloguing and photometry 

The NB 1060 and Y images are registered using stellar point sources 
as a reference. Similarly, the images are PSF matched by smoothing 
the Y image with a gaussian kernel until the stellar FWHM values 
are equal to those measured in the NB 1060 image. Source Extractor 
(Bertin & Amouts, 1996) is then run on the NB 1060 image to pro- 
duce a NB1060 selected source catalogue. Extraction parameters 
are given in Table |2] We find slight deviations in the background 
of the order of a few percent in the NB1060 detector 2 image. This 
appears to be due to radio active decay tracks created near to the de- 
tector. For this reason, for this detector only, the background is de- 
termined by filtering over smaller areas to produce a more accurate 
local background for each object. The altered parameters are given 
in the final column of Table |2] Values were chosen to minimise the 
fluctuations in the background measured by Source Extractor (as- 
sessed by outputting a 'check image' of type 'BACKGROUND') 

For each object in the catalogue, two flux measurements are 
made: K -A'SlOeO colour and total NB1060 flux. Circular aper- 
tures are used throughout. The size of each aperture is tuned to the 



size of the object as described in Labbe et al. (2003fl The Labbe et 
al. scheme recommends different sized apertures for making colour 
and total flux measurements, including modifications for blended, 
extended and particularly compact objects. Colour measurements 
are made in circular apertures of diameter DcoL=2(A/5o/;r)^/^, 
where Aj^o is the measured isophotal area within the detection 
isophote. Similarly, total flux measurements are made in apertures 
of diameter T)joT=^(Ai^i^oN / where A^ron is the Kron area; 
the area of the ellipse defined by the Kron radius (see ,Kroni,1980l) 
(in Source Extractor this is the area of the 'AUTO' aperture). 

For colour measurements, the object is detected in the narrow 
band and then equal sized circular apertures are placed in the same 
position on both the NB1060 and Y images, in an analogous way to 
using Source Extractor in dual image mode. Colour measurements 
therefore have the same spatial origin. 

Objects within ~ 10" of the edges of the stacks are removed 
due to the lower exposure time in these regions, leaving a survey 
area of 46.24 square arcmin. Detections with a S/N < 3.0 in the 
narrow band are additionally removed from the catalogue. 

There are 2150 objects in the full catalogue. We find the 5cT 
NB1060 flux limit in a circular aperture, 10 pixels (1.06") in diam- 
eter is 1.5 X 10^'^erg cm^^s^' (m^gio60=24.55). To estimate the 
point source completion of the catalogue, we insert artificial point 
sources into the images using the IRAF program MKOBJECT and 
extract them using Source Extractor as described above. We find the 
90% completeness limit of point sources in the NB 1060 catalogue 
ismA,Bio60=24.40(1.73 x 10^'^erg cm^^s"'.) 



3 EMISSION LINE GALAXY SELECTION 

We expect emission line galaxies (ELGs) to have excess NB1060 
flux compared to the Y band continuum. ELGs are therefore se- 
lected based on a clear flux excess in the narrow band, ie, Y — 
NB1060 > 0. An objec t is selected based o n two criteria: 

The parameter Z teunkeretalll 19951) is used to characterise 
the significance of the NB1060 excess compared to a flat spec- 
trum, taking into account the noise properties of the images. For 
this work, the appropriate selection curve is given by 



my - rrif^B '- 



-2.51ogi, 



_Y_lQ-QACiQ.O-mnB) 



NB ' 



(1) 



where (Jnb and Oy are the noise in the NB1060 and Y images re- 
spectively and 30.0 is the scaled zero-point as described in Section 
[23] 

We use a co lour significance of Z=3 to select the ELG can- 
didates teunker e t al. ( 1995)). This assures that the fraction of or- 
dinary, non emission line galaxies scattered into the sample due to 
noise is very low, ~ 1 in 1000 objects. 

Secondly an observed equivalent width (EWq/,,) criterion of 
EW > 50A is imposed, equivalent to a colour cut of 



y-A'sioeo > 0.37. 



(2) 



This corresponds to a rest-frame equivalent width of 17. SA for 
[Oil] atz=1.85. 

Figure|2]shows the colour magnitude selection diagram along 
with the S=3.0 selection line (curved) and the equivalent width cri- 
terion (solid horizontal line). Candidate ELGs are highlighted in 

' We note that if Source Extractor's elliptical 'AUTO' aperture is used in 
place of the TOT aperture presented here, this changes the final SFR mea- 
surements by < 3%. 
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Figure 2. Colour-magnitude candidate selection diagram. The dashed red 
line shows the E=3 selection criterion and the solid red line indicates the 
50A observed (17. 5A rest frame) equivalent width cutoff. ELG candidates 
are highlighted in black. 



Table 3. Emission Lines Possibly Detected in the NB1060 filter 



Emission 
Line 


(Mm) 


z 


dz 


(Mpc"3) 


log(LL„„)^ 

erg s"' 


M,,, yr"' 


Ha 


0.6563 


0.62 


0.016 


950 


40.38 


0.19 


[OIII] 


0.5007 


1.12 


0.021 


2353 


41.02 




[OIII] 


0.4959 


1.14 


0.021 


2410 


41.04 




H/3 


0.4861 


1.18 


0.021 


2528 


41.08 




[Oil] 


0.3727 


1.85 


0.028 


4138 


41.56 


6.34 


Lya 


0.1216 


7.73 


0.085 


7977 


43.04 





" Co-moving volume defined by the NB 1060 filter width given in Table[T] 
^ Luminosity corresponding to a NB 1060 line flux of log(//,„j) = — 16.82 erg 
s"' cm the 5fT NB1060 limit of the data (1" aperture) 

Limiting SFR corresponding to the limiting luminosity in column 6, assum- 
ing SFR=7.9xlO-*2 Lna (Kennicutt 1998) and [OII]/Ha=0.45 (see Section 



ED) 



black. Using these criteria, 58 objects are selected as ELG candi- 
dates. 



3.1 Selection of [Oil] emitters at z=1.85 

The 58 NB selected ELGs are expected to be comprised of 
[OII]A3727, Ha, [OIII]114959,5007 and H/3 emitters, along with 
a small fraction of redshift interlopers. Possible emission lines and 
their corresponding redshifts and volumes are given in Table[3] 

We split the sample into composite line samples in two ways: 
(1) We use the spectroscopic and photometric redshifts reported in 
the GOODS MUSIC catalogue, version 2 (Santini et al., 2009). (2) 
We use galaxy evolution tracks in the colour-colour diagram z — 7, 
y — ; to split the ELG sample into low and high redshift populations 
according to colour. 

Photometric redshifts in the MUSIC catalogue are based on 



photometry in 15 bands spanning the optical to the far infrared. 
We find good agreement between the 1697 spectroscopic red- 
shift measurements in the catalogue and their equivalent photo- 
metric redshift estimates. Splitting the MUSIC measurements into 
low (Zphot < 1-5) and high (Zphot > 1-5) redshift groups, we find 
the median offset of the low redshift group is {zpiwi — Zspec) = 
—0.004 with a dispersion of Omad = 0.072 (where (Jmad is the 
standard deviation associated with the median absolute deviation 
of {Zphot ^ Zspec}) and equivalently for the high redshift group, 
{zphot — Zspec) = +0.05 with a dispersion of Gmad = 0. 18. 

Our catalogue is matched to GOODS MUSIC using a search 
radius of (0.5")- All but two candidates have MUSIC matches. 
Figure [5] shows a histogram of the photometric redshift measure- 
ments (ZPHOT parameter in MUSIC) of candidate ELG galaxies. 
The histogram clearly shows 3 peaks corresponding to [Oil] A 3727 
at Zphoi ~ 2, Ha at Zphor ~ 0.6, and a merged peak at Zphor ~ 1 
comprised of [OIII]AA 4959,5007 and Hj3 emitters. We discard 
5 galaxies with Zphot < 0.3 or Zphot > 2.7 as redshift interlopers 
(~8% of the sample), noting that one with a photomet ric red- 
shift of 6.88 is likely to be a T-dwarf jEvles et al.ll2007h . Gaus- 
sian curves are fitted individually to each population. Fit parame- 
ters are Zphg, pll] ~ 2-05 ± 0.27 for the [Oil] candidate population 
and Zphot.Ha = 0-62 ±0.032 for the Ha population. The broader 
spread of the [Oil] peak compared to the Ha reflects the reduced 
accuracy of the photometric redshifts at high redshift. 

For the second diagnostic, we calculate galaxy evolution 
tracks in z — 7, V — i colour space us ing the photometric redshift 
code detailed in iBanerii et al.l ( |20I(]|) . Galaxy types E, Sbc and 
Scd are modelled using the average observed spectra of Coleman, 
Wu & Weedman ( 1980). In add i tion w e use the observed starburst 
model (SB2) from lKinnev et alJ h996) and a synthetic spectrum of 
a galaxy with a n instantaneous 50Myr star burst , generated using 
the PEGASE code jFioc & Rocca- Volmerange|[l997h . 

Each model is redshifted in steps of 0.01 producing tracks 
from z=0.6 (the lowest redshift. Ha emitters) to z=1.85 (the highest 
redshift [Oil] A 3727 emitters). 

In Figure |4l the evolutionary tracks are indicated by black 
lines, where (from top to bottom) the elliptical track is shown in 
the dotted line, type Sbc in dash, Scd in dash-dot, the 50Myr star- 
burst model in long-dash and the starburst (SB2) model in the solid 
line. Joining up the points on the five tracks at the redshift of each 
line emitter in the filter delineates where each population of line 
emitters is expected to lie in colour-colour space. The blue line 
joins points at z=1.85, indicating where [OII]A3727 emitters are 
expected to lie. Similarly, the red line joins points at z=0.63 corre- 
sponding to Ha and lines joining points at z=1.12,1.14 and 1.18 are 
shown in green corresponding to [OIII]AA4959,5007 and H/3 (top 
to bottom respectively). 

We use z — 7, V — i photometry from the MUSIC catalogue 
to over-plot the ELGs. The ELGs in the sample, indicated by cir- 
cles, are blue in z — 7, V — consistent with them being star form- 
ing galaxies. The candidates clearly split into two regions in colour 
space, the lower stream attributed to [Oil] objects and the upper to 
a combination of Ha, H/3 and [OIII] emitters. 

Final classifications, taking both colour and photometric red- 
shift into account, are colour coded on the Figure with objects clas- 
sified as [Oil] shown in blue (26 objects), Ha in red (14 objects) 
and H/3 or [OIII] objects in green (13 objects). One object has a 
photometric redshift placing it as an [OIII] or H/3 emitter, yet ap- 
pears in the [Oil] stream in Figure |4] Upon visual inspection, this 
object appears to be blended with a neighbouring object in the im- 
age which may be affecting the measurements and we therefore 
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Figure 3. MUSIC catalogue photometric redshifts for the sample of ELGs, 
showing three peaks corresponding to [Oil] emitters (blue). Ha emitters 
(red) and a combination of [OIII] and H/3 emitters (green). Interloping 
galaxies are shown in black. 

discard it from the sample. Ten objects have good quality spectro- 
scopic redshifts and these are indicated in the Figure with asterisks. 

Over all, there is good agreement between the track, photo- 
metric redshift and spectroscopic classifications. When the Zphoi 
classification is taken into consideration, objects in Figure|4]clearly 
split into three classes, with [OIII] and Hj8 emitters lying bluer than 
Ha in y — i. 

Four objects are left unclassified: two with no MUSIC match 
and a further two with partial entries in MUSIC, not including JHK 
or photometric redshift measurements. These are classified using 
B — i, V —z colours in the same manner as in Figure |4] For the two 
objects with no MUSIC data, colours are determined from ACS 
cutouts generated by the MAST cutout tool and the documented 
zero-points for the ACS images. On the basis of the positions of 
these objects in B — V — z space, three were classified as [Oil] 
emitters and the remaining object was classified as either an [OIII] 
or Hj3 emitter. 

We checked the sample was uncontaminated by stars by vi- 
sually inspecting each object and additionally ensuring that none 
of the objects were flagged as stars in the MUSIC catalogue. None 
of the objects were flagged as AGN in MUSIC and none of our 
objects matched the Chandra Deep Field-South: 2Ms Source Cata- 
logue (Luo et al. 2008), suggesting that the sample is additionally 
uncontaminated by AGN. 

In summary, of 58 initial ELG candidates, 53 were identified 
as genuine ELGs. Our final [Oil] emitter sample contains 26 ob- 
jects. 



4 THE OBSERVED LUMINOSITY FUNCTION 

The observed luminosity function is calculated in 4 stages as fol- 
lows: In Section [42] the luminosity function is computed under 
the simplifying assumption that the NB1060 filter is a perfect top 



Figure 4. Splitting ELGs into low and high redshift groups according to 
position in z — J,V — i space. Black lines show galaxy evolution tracks of 
different galaxy models redshifted from z=0.6-z=1.85. Top to bottom these 
are: E, She, Scd, 50Myr instantaneous starburst and Star Burst model SB2 
(See text for details). The blue line joins up points on the tracks coiTespond- 
ing to models at z=1.85, indicating where [Oil] emitters are expected to lie 
in the diagram. Similarly, the red line joins points at z=0.6, showing where 
Ha emitters are expected to lie and the green lines join points at z=1.12,1.14 
and 1.18 corresponding to [Oni]AA4959,50()7 and H/3. Candidate emitters 
are colour coded according to their photometric redshifts as in Figure|3]with 
[Oil] candidates in blue, [OIII] and H/3 in green and Ha in red. Asterisks 
indicate spectroscopically confirmed candidates (Santini et al. 2009). 



hat function with a width equal to the FWHM of the real filter. In 
Section l431 we investigate the completeness limit of the survey. In 
Section |4!4l we fit a Schechter function to the complete region of 
the luminosity function and finally in Section |43] simulations are 
undertaken to scale the luminosity function to take into account 
the real filter shape and produce the corrected observed luminosity 
function. 



4.1 Line fluxes and [OH] luminosities. 

The narrow band contains both line and continuum emission. We 
correct for continuum emission using the Y flux density. Given that 
the NB1060 and Y bands overlap, the appropriate equation for the 
line flux is given by: 

r fNBW60 - £fY 

.// = : (i) 

1 -e 

where f^BlQeQ and fy are the total fluxes in the narrow and broad 
band filters and e is the ratio of the widths of the narrow and broad- 
band filters (in this survey e=0. 102). 

Assuming all objects in the [Oil] sample lie at the centre of 
the filter at z=1.85 and the luminosity distance, d/, = 4.36 x 10^^ 
cm, the observed [Oil] line luminosities are: 

L[01lUs = ^Ttdlf,. (4) 
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Table 4. Properties of the [Oil] Sample. 
R.A^ DEC" NB^ a^if^ 



ay 



tog(/Lll 



iog(L|o;/|„fc^ 



a co-moving volume of Vc=4138 Mpc . The fixed- volume lumi- 
nosity function is tabulated in Table|5] 



(J2000) 


(J2000) 






















53 


.15124 


-27.79239 


23, 


.02 


0.05 


24, 


.39 


0.08 


-16 


.32 


42 


.06 


53 


.15497 


-27.79037 


23, 


.52 


0.08 


24, 


.35 


0.09 


-16 


.67 


41 


.71 


53 


.16162 


-27.78743 


22, 


.58 


0.04 


23, 


.26 


0.04 


-16 


.36 


42 


.02 


53 


.15561 
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" Measurements in degrees 
^ NB 1060 and Y magnitudes are on the AB scale 
Photometric En'or 



Line Flux (erg cm ' 



Observed Line Luminosity (erg s ) 



NB 1060 and Y magnitudes, along with the line fluxes and line 
luminosities for the sample of 26 [Oil] emitters are given in Table 
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4.2 Fixed-volume luminosity function 

The luminosity function is initially calculated by assuming the 
NB1060 filter is a perfect top hat function. Under this approxima- 
tion, all objects are visible through the FWHIVI of the filter, lead- 
ing to a constant survey volume. Additionally, we assume the filter 
is sufficiently narrow such that it has uniform sensitivity to line 
strengths throughout the full filter width (ie. moving an object of 
fixed intrinsic line luminosity across the redshift range defined by 
the filter will not significantly alter its observed line flux.) 

For clarity, the resulting binned luminosity function is referred 
to as the 'fixed-volume' luminosity function, given by the relation 



^i{logL{[01I])) - 



1 



A(logL([0//l)) Vc 



(5) 



where Vc is the fixed co-moving volume probed by the filter, 
A(log L([OII])) is the bin width and N,- is the number of galaxies 
with [Oil] luminosity in the range log L([OII])± 0.5Alog L([OII]). 
When assuming a top hat filter function, Vc is fixed to the volume 
defined by the FWHM of the NB1060 filter (0.0104 /xm.) With a 
survey area of 46.24 square arc minutes, at z=1.85 the survey covers 



4.3 Sample completeness 

4.3.1 Detection limit 

By using the NB1060 image as the detection image (Source Ex- 
tractor in Dual Image IVlode), it is not necessary for an object to be 
detected in the Y band for it to be included in the sample. In princi- 
ple therefore, the survey is sensitive to objects with infinitely high 
equivalent width lines. 

Given that objects are detected solely in the naiTow band, the 
limiting line flux is closely linked to the NB1060 detection limit. 
The detection limit and detection efficiency are explored by intro- 
ducing synthetic populations of objects into the images using the 
IRAF package MKOBJECT. These are then recovered using the 
detection techniques described in Section l23l The detection com- 
pleteness is the percentage of input objects that are successfully 
extracted from the images. 

As expected for high redshift populations, the real [Oil] emit- 
ters are not well resolved in the HAWK-I images. The average mor- 
phology of the [Oil] emitters (measured from the HAWK-I stacked 
image) is approximated by an exponential disk with scale length 
r5=2.5". Due to the high redshift of the objects and the low reso- 
lution, we find the average profile to be representative of the [Oil] 
population as a whole. For this reason, the simulations were limited 
to a single input profile. 

Batches of 100 objects were introduced into the NB1060 im- 
age with continuum magnitudes 18-27 in steps of 0.2 mag. Source 
Extractor was run on the images using the detection parameters 
given in Table |2] and the detection rate was measured. The results 
of the simulations can be seen in Figure|5] which shows how the de- 
tection efficiency of the [Oil] emitter profile varies with input nar- 
row band magnitude. The Figure indicates that for the [Oil] emit- 
ter profile, the survey is 90% complete to magnitude 24.2 (cf. 90% 
completeness to magnitude 24.4 for point sources (Section [2.3t ). 



4.3.2 Line flux completion 

Whilst sensitive to high equivalent widths, the equivalent width 
threshold applied in the selection procedure in Section |3]results in 
objects with EWq^, < 50A (equivalent to EW^,,; = 17. 5 A) being 
omitted from the survey. 

Figure|6]shows the selection diagram. Objects above the Z = 3 
selection line are highlighted in bold. ELGs that are potentially 
missed by the EWq/^, > 50A threshold lie between the curved and 
solid red lines. Of the objects falling in this category, we find most 
of them are either stars (highlighted in orange asterisks) or spec- 
troscopically confirmed redshift interlopers (red crosses). Such a 
high proportion of confirmed redshift interlopers indicates that the 
equivalent width threshold chosen in this study is largely robust. 

Of the remaining objects, only five have MUSIC photometric 
redshifts that could lead to them being classified as [Oil] emitters, 
these are highlighted with blue diamonds. We assess the signifi- 
cance of this group of omitted objects by over-plotting lines of 
constant [Oil] line luminosity, representing the edges of the bins 
of the fixed- volume luminosity function in Table |5] For a bin to 
be complete, all objects between the lines of constant line lumi- 
nosity corresponding to the edges of the bin must be selected. For 
example, for the bin centred on log(L) = 41.74 erg s^' to be com- 



[Oil] emitters in the GOODS field at z~ 1.85; a homogeneous measure of evolving star formation. 
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nput NB1060 (mag) NB1060 (AB) 



Figure 5. Detection completeness as a function of magnitude for synthetic 
galaxies modelled with the average [Oil] emitter profile, an exponential disk 
with rj=2.5". 

plete, all objects between the lines log(Lj(5//]) = 41.64 erg s^' and 
log(LjQ//]) = 41.84 erg s^' must be selected. 

It can be seen that three of the five potential [Oil] candidates 
below the equivalent width threshold have colours placing them 
in the two lowest luminosity bins in the fixed-volume luminosity 
function. These bins both lie fainter than the 90% detection limit of 
i^NBlOdQ = 24.2 defined in Section l4.3.1l The two lowest luminosity 
bins in the study are therefore incomplete due to a combination of 
detection and selection incompleteness. 

The small number of possible [Oil] emitters falling in brighter 
bins leads us to conclude that we do not lose a significant fraction 
of emitters using this EW criterion. 

Furthermore, completeness simulations (using the [Oil] pro- 
file described in Section 14.3. It show a recovery rate greater than 
90% for [Oil] emitters with Y band continuum magnitudes fainter 
than magnitude 23 and lines brighter than log(L/,„(.)=41.74. This 
recovery rate rapidly drops off for fainter lines. These results, to- 
gether with those presented in Figure |6] lead us to conclude that an 
appropriate completeness limit for the study is the minimum lumi- 
nosity, log(L„„„) = 41.74 erg s^' . 

4.4 Schechter fit 

The fixed-volu me luminosity function is fit with a Schechter func- 
tion (Schechter 1976) of the form 

To compute an appropriate Schechter function for the bright 
end of the luminosity function, a range of faint end slopes are as- 
sumed. 

For comparability, we vary a over the same range as that as- 
sumed by Zhu et al. (2009) in their [Oil] survey at 0.75 < z < 1 .45: 
a = — 1.3±0.2. For comparison, Takahashi et al. (2007) measured 



Figure 6. Same as Figure|2] overlaid with lines of constant line luminosity, 
labelled logarithmically. The luminosities chosen represent the edges of the 
bins used to create the fixed- volume luminosity function in Table |5] A bin 
is complete if all the objects between the bin edges are included in the sam- 
ple. Candidate [Oil] emitters are highlighted in blue, stars are indicated by 
orange asterisks and objects with redshifts placing them out of range of the 
filter range are shown in red. 

the faint end slope of the [Oil] luminosity function at z= 1.2 in two 
fields finding a = -l-41^o [5 and a = -1.38+0 3°. 

In Section |5] we look at how the luminosity function evolves 
from 2=1.2-1.85. We emphasise that the range of a assumed here 
has a minimal effect as this analysis is restricted to the bright por- 
tion of the luminosity function that is well fit, regardless of the 
assumed faint end slope. 

We fit the [Oil] luminosity function at L > L,„,„ with a 
Schechter function, assuming a = —1.3 ±0.2, using the maxi- 
mum likelihood parametric fit method (Sandage, Tammann Sl Yahil 
1979). The resulting best fit parameters are log(L*) = 42.00 ± 
0.06 erg s"' and log(0*) = -2.21 ±0.09 IVIpc-^ We note that in- 
creasing Linin produces fits within these eiTors but lowering L,„,„ 
quickly departs from them, indicating that log(L„„„) = 41.74 erg 
s^' is an appropriate completeness limit for the survey. 

4.5 Filter correction 

The fitted fixed-volume luminosity function is scaled to take into 
account the effects of filter shape. There are two considerations: 

(i) Flux loss due to filter transmission: ELGs with [Oil] lines 
that fall in the filter wings have lower observed line fluxes due to 
the poor transmission. When computing the luminosity function, 
this results in a proportion of line objects systematically moving 
into fainter bins. 

(ii) Variation of volume with line strength: Whilst a relatively 
faint emission line in the filter wings may fall below the detection 
threshold, a bright line will still be detectable. Brighter lines are 
therefore detectable over a wider filter width and correspondingly 
over a wider redshift range and a larger volume. 

Simulations were run to quantify these effects. We considered 
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Table 5. The Fixed- Volume Luminosity Function 
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Observed Luminosity Function 


log L([01I]) 


log(j) 
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41.34 


-2.22 
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-2.14 
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-2.14 


6 


42.14 


-2.44 


3 


42.34 


-2.92 


1 



a volume large enough to fully encompass the filter volume, includ- 
ing the filter wings. The volume was populated with objects with 
density and flux distributions according to a trial input luminosity 
function. We assumed the objects were homogeneously distributed 
with respect to redshift. The filter profile was then used to recover 
the observed simulated object luminosities. Different values of in- 
put (j)* and L* were iterated through until the output of the simula- 
tion matched the observed [Oil] luminosity distribution. 

The simulations indicate that the intrinsic observed luminos- 
ity function, assuming a = —1.3 ±0.2 is best fit with Schechter pa- 
rameters of log(L*) = 42.05 ±0.06 erg s"' and log(0*) = -2.23± 
0.09 Mpc"^ 

Figure |7] shows the observed [Oil] luminosity function found 
in this study alongside the observed luminosity functions of equiv- 
alent [Oil] studies at z=1.2. Black triangles show our binned lu- 
minosity function, calculated assuming a constant filter volume. 
The luminosity above which the survey is estimated to be complete 
(Lii„) is indicated by the vertical dotted line and the two bins that 
are significantly incomplete are highlighted as lower limits. The 
best fitting Schechter function to the complete portion of the lumi- 
nosity function is indicated with the black dashed line. Finally, the 
solid black line is the filter-corrected luminosity function. Over- 
plotted are comparable results from other [Oil] surveys at z=1.2. 
In red is the Schechter fit of Ly et al. (2007), green points show 
the binned luminosity functions of Takahashi et al. (2007) in the 
COSMOS field (diamonds) and SDF (asterisks) and the binned lu- 
minosity function of Zhu et al. (2009) is shown in blue. Note that no 
correction for dust obscuration has been made to any of the lumi- 
nosity functions presented in the figure. An appropriate obscuration 
correction for our sample is calculated in Section |6j2] 



5 EVOLUTION OF THE [Oil] LUMINOSITY FUNCTION. 

Figure |7] suggests the observed [Oil] luminosity function evolves 
between redshift 1.2 and 1.85. To quantify this evolution, we com- 
pute the number density of [Oil] emitters at a range of redshifts. 
We concentrate on the integrated number density of objects in 
the luminosity range in which our luminosity function is robust: 
log(V/]) > 41.74 ergs- 1. 

We find the number density of objects in the luminosity 
range log(Lj(5//]) > 41.74 erg s^' at z=1.85 in this survey is 
log(Plog(L|o„|)>4l.74) = -2.45 ±0.14 Mpc-3. Figure [8] shows this 
result alongside the equivalent number densities of Zhu et al. (2009) 
(z=0.84, 1.00, 1.19 and 1.35) in blue squares, Ly et al. (2007) 
(z=0.91 and 1.18) in diamonds and Takahashi et al. (2007) (z=1.2) 
in red crosses. Our measurement is indicated by the red triangle. 

We fit the points (all but the two outliers) with a line. 
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Figure 7. Comparison of our observed luminosity function at z=1.85 with 
other [Oil] luminosity functions from the literature at 2=1.2. Black triangles 
show the binned fixed-volume luminosity function from this work, with the 
best fitting Schechter function in the black dashed line. The final observed 
fit, coiTected for filter shape is shown in the sohd black line. The HAWK- 
1 hmiting luminosity is indicated by the vertical dotted line. The observed 
luminosity functions reported by Takahashi et al. (2007) are shown in green, 
that of Zhu et al.(2009) in blue and the Schechter fit of Ly et al. (2007) is 
shown in red. Note that none of the luminosity functions shown here have 
been corrected for obscuration. 
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Figure 8. Evolution of the observed [Oil] luminosity function: evolution 
of the total number of bright objects per Mpc^^ in the luminosity range 
log(Ljo;;|„j,) > 41.74. The dashed lines indicate the best and ±a fits given 
by the equation log(<I>[log(L[o,/]„j„) > 41.74,z]) =mz + c. 



log(*[log(L[o„]) > 41.74, z]) = mz + c, finding in = 0.69±0.089 
and c = —3.72 ±0.11, this is over-plotted in green on the Figure, 
along with the ±(T fits. We find that the space density (p) of bright 
objects at z=1.85 in the luminosity range log(L|Q//|) > 41.74 erg 
s^' is a factor 2 greater than the observed space density of [Oil] 



[Oil] emitters in the GOODS field at z~ 1.85; a homogeneous measure of evolving star formation. 
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emitters reported at z~1.4; a comparable increase to that reported 
between z=1.0 and z=1.4. 



6 THE STAR FORMATION RATE AT Z=1.85 

Two assumptions have to be made to convert the [Oil] luminosity 
function found in this work into a star formation rate density: (i) 
The [Oil] luminosities have to be corrected for internal extinction 
and (ii) a calibration has to be assumed between [Oil] luminosity 
and SFR. 

In Section [6T| we explore the use of a "common" obscura- 
tion correction to make the results comparable to the compilation 
of SFR measurements made by Hopkins (2004). In Section[6j2]we 
extend the analysis and independently estimate the extinction cor- 
rection using 24/im flux and discuss the appropriate [OII]/Ha ratio. 
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6.1 Common obscuration correction 

iHopkinj ( l2004h provides a compilation of SFR measurements 
made using a range of SFR indicators, corrected according to a 
common framework. In order to compare the measurement made 
here to the Hopkins compilation, we convert our luminosity den- 
sity to SFR according to this common framework. 

The first step towards computing the SFR is to calculate the 
total [Oil] luminosity density, Jif[OII], at z=1.85 by integrating the 
Schechter function: 



(7) 



^[OII]= f <^(L)LdL = ^*L*T(a + 2) 
./() 

where F represents the gamma function. 

Hopkins converts the [Oil] luminosity density to Ha assum- 
ing /.2//a,oij=0.45. The inferred Ha luminosity is then cor- 
rected for internal extinction, assuming A/fc(=1.0 mag. Note that 
A//a=1.0 mag corresponds to Aj(j;/] = 1.86 using the O'Donnell 
(1994) galactic obscuration curve (R,,=3.1). Finally 
converted in to a SFRD using the calibration of Kennicutt (1998): 



: 7.9 X 10-42 _^Ha,coiT Meyr-^ 



(8) 



The [Oil] Luminosity function derived in this work, corrected 
for filter effects (Schechter Parameters: a = — 1.3±0.2, log(L*) = 
42.05 ±0.06 ergs"' andlog((|)*) = -2.23 ±0.09 Mpc"^), implies 
an [Oil] luminosity density of log(^[0//]) = 39.93 ±0.08 erg s"' . 
Accounting for A//a=lmag, this implies a SFRD of = 0.38 ± 
0.06 Mq yr^'Mpc^^ 

This result is plotted in red in Figure|9] alongside other SFRD 
measurements derived from [Oil] and corrected using the "com- 
mon" scheme at lower redshift. 

Whilst the application of a common obscuration correction 
is helpful for comparison with published literature (see Hopkins, 
2004), it is an oversimplification. A more rigorous analysis can be 
carried out by measuring the obscuration directly. This is calculated 
and discussed in detail in Section l631 The final SFRD calculated in 
Section |6^ is shown for completeness in orange in Figure|9] 



6.2 Independent measure of A.^on] [OII]/Ha 

To independently estimate A.^oil] f*"" our sample, we look to the far 
infrared. Kennicutt et al. (2009, hereafter K09) developed a set of 
empirical calibrations converting [Oil] line luminosities into SFRs 
using 8, 24/xm or total infra-red (TIR) fluxes as a tracer for the dust 
emission. They consider two samples of z=0 galaxies; the SINGS 



Figure 9. A compilation of SFR measurements using the [Oil] line, as- 
suming a constant obscuration correction, AHa=l mag. [Oil] measurements 
from the compilation of Hopkins (2004) ai'e shown in grey, overlaid with the 
measurements of Hippelein et al. (2003) (purple), Takahashi et al. (2007) 
(green) and Zhu et al. (2009) (blue). The red triangle shows the result from 
this work corrected using the same "common" correction as applied to the 
lower redshift points (AHa=l mag). The orange triangle shows the same 
result, this time corrected using Ajq,;] = 0.98 mag (Ana = 0.52), the mean 
obscuration of the z=1.85 sample inferred from 24/im fluxes (see Section 

E3. 



sample of 75 local galaxies with distances less than 30 Mpc, pre- 
sented in Kennicutt et al. (2003) and 417 galaxies from the survey 
of integrated spectrophotometry described in Moustakas & Ken- 
nicutt (2006, hereafter IVIK06). The samples were selected to be 
representative of the wide range of morphologies, luminosities and 
dust opacities seen in present day galaxies. The combined sample 
includes objects ranging from dwarf irregulars to giant spirals and 
IR-luminous galaxies. Full details can be found in the respective 
survey papers. For these local galaxy samples, they find the cor- 
rected [Oil] luminosity density is given by 

L[OII\c„rr = L[OIlUs + 0.016L(8/im) (9) 
implying 

0.016L(8/xm),.e 



'[0//1 



= 2.5 log 



L[01I] 



obs 



(10) 



For galaxies at the redshift of this survey, 8jUm flux is red- 
shifted to ~ 24/im. Using the 24/im fluxes from the IVIUSIC cat- 
alogue, we calculate values of A.^oil\ sample. MUSIC con- 
tains measurements for 24 of the 26 [Oil] emitters in the survey, 16 
of which are upper limits. 

To compare our results with galaxies at z=0, we derive the 
[Oil] obscuration values for the K09 and MK06 samples using 
the stellar-absorption-corrected Ha/Hj8 ratios quoted in K09 and 
MK06. The ratios are converted to [Oil] obscurations using the 
same assumptions as detailed in K09; namely we assume an intrin- 
sic Ha/H/3 ratio for Case B recombination of /(Ha)//(H/3)=2.86 
(electron temperature T^ = 10,000 K and density Ng = 100 cm^-^). 
The observed reddenings are then converted to [Oil] attenuation 
values (via Ha) using the O'Donnell (1994) extinction law, assum- 
ing Rv=3.1. 

Figure [Tol a) shows the distribution of Ajq//] with respect to 
the uncorrected [Oil] luminosity. The z=0 SINGS and MK06 sam- 
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Figure 10. (a) [Oil] extinction witli respect to observed [Oil] luminosity for galaxy samples at 2=0 (MK06 and SINGS samples shown in diamonds and crosses 
respectively) and z=1.85 (HAWK-I shown with blue triangles). 16 of the 24 A|o;;| at z=1.85 were calculated from upper limit measurements of the 24/im 
flux (indicated by the arrows). The eight measurements made from positive detections at 24/im are highlighted with squares for clarity. Red and Blue solid 
lines indicate the mean A|o;;| values of the combined MK06 and SINGS sample and the Hawk-I sample respectively. Dotted lines indicate the en'or on the 
mean in each case, (b) [Oil] extinction with respect to obscuration corrected [Oil] luminosity, where a correction has been made according to the measured 
obscuration values. The green line shows the maximum obscuration with respect to intrinsic luminosity that could be selected in our sample. 



pies are shown in crosses and diamonds respectively. The z=1.85 
results from this work are over-plotted in blue triangles. For the 
z=1.85 objects in our sample, values of A|(j//] derived from upper 
limits are indicated with arrows and the eight values of A|o//] de- 
rived from positive detections at 24/im are highlighted with squares 
for clarity. 

Although the z = measurements of A.\on\ have a large dis- 
persion, we find no evidence for a systematic variation of f^^oil] 
with observed [Oil] luminosity. The mean of measurements with 
log(LjQ//]„ii,j) > 40.0 is AjQ,/] = 1 .59 ± 0.07 mag with a dispersion 
of 0.8 mag. Bifurcating the sample into 40.0 < log(L[o//]ofcj) < 41 .0 
and 41.0 < log(i[o//]oi.s) < 42.0 yields mean values of Aj^//] = 
1.42 ±0.1 and Apjj^ = 1.80 ±0.10 (dispersions of 0.8 and 0.7 
mag) for the fainter and brighter samples respectively, consistent 
with no luminosity dependence. The mean values of each sample 
are indicated by the red solid lines in Figure [lO]and the red dotted 
lines indicate the error on the mean in each case. 

Our galaxy sample at z=1.85 however has systematically 
lower levels of Ajq//] than the local universe galaxies, with a mean 
^[Oll] ~ 0-98 ±0.11 mag and dispersion of 0.6 mag. This corre- 
sponds to A/fa = 0.52 (O'Donnell (1994) extinction curve assum- 
ing Ry =3.1). 

Selection effects may bias this result. For any given intrin- 
sic luminosity, the survey is biased towards selecting objects with 
low obscuration: objects with high obscuration values may have 
(observed) [Oil] fluxes that fall below the selection threshold of 
the survey. To investigate the effect this would have on the mea- 
sured mean Ajq/;j value, we re-plotted Figure [ToT a), coiTecting the 



observed luminosities for extinction using the measured obscura- 
tion corrections, this is shown in panel (b) of Figure [TO] The green 
line shows the maximum detectable value of A[OII] as a func- 
tion of corrected-luminosity. The mean obscuration of the SINGS 
and MKQ6 galaxies with observed [Oil] fluxes brighter than the 
Hawk-I selection threshold is A[OII]= 1.68, whereas the mean 
of all objects with \ogL[OII]co,r > 41.5 = 1.90. Therefore, ap- 
plying the Hawk-I selection bias to the z=0 sample would result 
in a ~ 12% drop in the mean A[OII] value. It is possible that 
this selection limit imposes more than the locally inferred 12% 
bias on the z=1.85 sample if a greater proportion of star forma- 
tion is dust enshrouded at high redshift. This cannot be ruled out 
on the basis of the present study, although it is worth noting that 
we found no evidence for a population of bright, dust enshrouded 
[Oil] emitters (Section [4.3. 2t . Adopting the mean value of Ajq;/] 
at z=1.85 (A[o„] = 0.98±0.11, Ana = 0.52) yields a SFRD of 
O.24±O.O6M0 yr"' IVIpc"^ This is plotted in orange in Figure |9] 
We have indicated this result as a lower limit to reflect the fact that a 
proportion of dusty [Oil] emitters may be omitted from the sample 
due to the selection criteria. 

In the above analysis we have assumed the widely used 
[OII]/Ha ratio of 0.45 (Kennicutt 1992, 1998). However, local 
surveys have suggested that the [OII]/Ho! ratio is luminosity- 
dependent (Jansen et al. 2001) as well being dependent on metal- 
licity and obscuration (Kewley et al. 2004). Hopkins et al. (2003) 
noted higher [OII]/Ho! ratios in higher equivalent width sys- 
tems. For a complete sample of 752 SDSS galaxies, they found 



[Oil] emitters in the GOODS field at z~ 1.85." a homogeneous measure of evolving star formation. 
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[OII]/Ha = 0.23 but that this rose to 0.46 if an EW limit of 
EW(Ha) > 70yl was imposed. 

Applying the Jansen et al. (2001) empirical relation between 
the the rest-frame absolute B band magnitudes and the OII/Ha to 
our galaxy sample, we find a mean absolute B band magnitude of 
Baftj=-20. 1 mag (/ = 24.5 mag), corresponding to [OII]/Ha = 0.48 
with an rms dispersion of 0.1. This is very close to the value 
(0.45) that we have assumed, and would result in a SFRD of 
O.23±O.O6M0 yr"' Mpc"^ 

The discussion above highlights that the "common" obscura- 
tion correction is an oversimplification of the problem. Upon mea- 
suring the obscuration of the objects we find a mean obscuration 
of Ajo//] = 0.98 mag (A//a = 0.52), rather than A.Ha = 1 as as- 
sumed in the common framework. However, given that we also 
find a proportion of dusty emitters are omitted from the sample, 
the SFRD measurement based on the measured obscuration correc- 
tion is likely to be a lower limit. It is reasonable to expect that the 
reality may be somewhere between the two measurements. Deeper 
data would enable this to be explored in more detail. 



7 SUMMARY AND CONCLUSIONS 

This study has used Science Verification Data from the ESO instru- 
ment HAWK-I to perform a high redshift survey for [Oil] emitters 
in the GOODS field. The [0II]13727 doublet is of particular inter- 
est to SFRD surveys in that it is visible to z~ 5 in the near-infrared, 
compared to the z~ 2.5 limit of Ha surveys. Recent advances in 
the calibration of the [OII]A3727 doublet (see eg. Kennicutt et al. 
2009) have made it possible to use [Oil] with greatly improved 
precision, facilitating homogeneous measurements of the SFRD of 
the universe to be measured from z=0 to z=5 for the first time. At 
z=1.85, this is the highest redshift [Oil] survey to date. We have 
identified 26 [Oil] emitters in a volume of 4138 Mpc^, to a 5tT flux 
limit of 1.5 X 10^ '^erg cm^^s^'. Our findings can be summarised 
as follows: 

• The observed [Oil] luminosity function at z=1.85 can be fit 
by a Schechter function with log(L*) = 42.05 ±0.06 erg s^' and 
log((^*) = -2.23±0.09 Mpc"3^ assuming a = -1.3±0.2; a rep- 
resentative range of high-z faint end slopes, including the value re- 
ported by Takahashi et al. (2007) at z=1.2. This range was also as- 
sumed by Zhu et al. (2009) in their survey of z=1.2 [Oil] emitters. 

• The space density (p) of bright (log(L[on]obs) > 41.74) [Oil] 
emitters at z=1.85 is log(p) = -2.45±0.14 Mpc"-*: a factor of 2 
greater than the observed space density of [Oil] emitters reported at 
z~1.4. This is a comparable increase to that reported between z=1.0 
and z=1.4 and we find that the [Oil] number density evolution of 
objects in the range L|on]obs > 41.74 can be fit by the function, 
log(<5[log(L|ou|) > 4I.74,z]) = mz + c, finding m = 0.69 ±0.089 
andc= -3.72 ±0.11. 

• We convert the [Oil] fluxes into a SFRD using the 
"common" extinction correction (A//a=I.O) and [OII]/Ho! ratio 
([OII]/Ha=0.45) employed by Hopkins (2004), finding a SFRD at 
z=1.85 of p* = 0.38 ± 0.06 yr^' . When compared to other re- 
ported values of the SFRD, calculated using the [Oil] emission line 
and the "common" conversion, our work suggests a three fold in- 
crease in the SFRD between z=1.4 and z=1.85. 

• We independently estimate Aj^//] for each object using rest 
frame 8/Jm flux (observed 24^m) and the empirical calibrations of 
Kennicutt et al. (2009). The results indicate that the [Oil] emitters 
we detect contain low levels of dust - the mean extinction of the 



sample being A.^oil\ ~ 0-98 ±0.11 with a dispersion of 0.6 (equiv- 
alent to A//a = 0.52). This is almost half the value measured at 
z=0. 

• Possible explanations of the low dust content are: (i) Selection 
Effects: The survey could be missing bright, dusty [Oil] emitters 
due to the equivalent width threshold applied to the sample. This 
is unlikely as we find no evidence to indicate the survey misses 
bright objects with low-equivalent widths (Section l4.3.2t . We do 
find however that selecting objects above a fixed absented [Oil] 
flux threshold may account for some of the difference between 
our A[OII] measurement at z=1.85 with respect to z=0. This is be- 
cause objects with high levels of obscuration fall below the ob- 
serx'ed [Oil] flux threshold leading to an underestimation of the 
true value. We estimate this bias results in an underestimation of 
the true mean by at least ~ 12%. However, without deeper data, it 
is unclear whether the bias alone could account for the full (~ 50%) 
reduction in dust obscuration which is seen between the z=1.85 and 
unbiased z=0 samples, (ii) Cosmic variance: the survey is relatively 
small (~ 4000 Mpc^"*) and therefore may not be representative of 
a typical region at z=1.85. (iii) Reddening Law: The K09 empirical 
calibrations used to derive A.^oil\ at z=L85 in this work were mea- 
sured at z=0. Different reddening laws may apply at z=1.85. (iv) 
Lower dust content at z= 1.85: the tight dispersion supports the idea 
that there is genuine, measurable difference between the amount of 
dust in high and low redshift [Oil] emitters. This would fit in with 
previous reports of low dust content in high-redshift galaxies (see 
e.g. Bunker et al. 2010, Ho et al. 2010 and references therein). 

• Incorporating the high-redshift value of A.^oil] ~ '^•^S into the 
SFR estimate yields a final SFRD of 0.24 ± O.O6M0 yr"' Mpc^^ 
This is a lower limit on the star formation density at z= 1.85, given 
that a proportion of dusty emitters are omitted from the sample due 
to the bias noted above. This is the first result tracing the SFRD to 
z=1.85 using [Oil]. It is in agreement with the UV SFRD measure- 
ments of Reddy et al. (2008) who found p» = 0.21 (with an error 
~ 20%) at z ~ 2. 
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